Conformal Invariance of Iso-height Lines in two-dimensional KPZ Surface 
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The statistics of the iso-height lines in (2+l)-dimensional Kardar-Parisi-Zhang (KPZ) model is 
shown to be conformal invariant and equivalent to those of self-avoiding random walks. This leads to 
a rich variety of new exact analytical results for the KPZ dynamics. We present direct evidence that 
the iso-height lines can be described by the family of conformal invariant curves called Schramm- 
Loewner evolution (or SLE K ) with diffusivity k = 8/3. It is shown that the absence of the non-linear 
term in the KPZ equation will change the diffusivity re from 8/3 to 4, indicating that the iso-height 
lines of the Edwards- Wilkinson (EW) surface are also conformally invariant, and belong to the 
universality class of the domain walls in the 0(2) spin model. 
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Recently, it was shown that the statistics of the zero- 
vorticity lines in inverse cascade of two dimensional (2D) 
Navier-Stokes turbulence is conformally invariant and be- 
longs to the percolation universality class The same 
issue has been studied for zero-temperature isolines in the 
inverse cascade of surface quasigeostrophic turbulence Q , 
domain walls of spin glasses [3( and the nodal lines of 
random wave functions 0]. Moreover, it has been shown 
recently that the statistics of the iso-height lines on the 
experimental WO3 grown surface is the same as domain 
walls statistics in the critical Ising model as well as those 
of Ballistic Deposition (BD) model 

Evidence of conformal invariance in the geometrical fea- 
tures of such complex nonlinear systems have been pro- 
vided, in the continuum limit, by stochastic (Schramm) 
Loewner evolution, i.e., SLE K , where k is the diffusivity 
[1, 0] • Schramm and Sheffield showed that the contour 
lines in a two-dimensional discrete Gaussian free field are 
statistically equivalent to SLE4 Q. Moreover, it is shown 
that the restriction property only applies in the case for 
k = 8/3 1. Since self-avoiding random walk (SAW) sat- 
isfies the restriction property, it is conjectured in the scal- 
ing limit to fall in the SLE class with k = 8/3 [TO'l The 
scaling limit of SAW in the half-plane has been proven 
to exist [ll[ but there is no general proof of its existence. 

In this Letter we investigate numerically the iso- 
height lines of the (2+ 1 )-dimensional Kardar-Parisi- 
Zhang (KPZ) model [12||, and study their possible con- 
formal invariance. It is shown that the KPZ's iso-height 
lines are equivalent to self-avoiding walks, and that the 
iso-height lines in the 2D-KPZ surface are SLEg,/ 3 curves. 
For the Edwards- Wilkinson (EW) interface (the KPZ 
model without the nonlinear term) the iso-height lines 
fall in the universality class of the interfaces in the 0(2) 
model, and can be described by SLE4. 

The KPZ equation is given by 

^|^--V^+^|V/ l | 2 +,(x,t). (I) 

The first term on the r.h.s describes relaxation of the in- 
terface caused by a surface tension v, and the nonlinear 



term is due to the lateral growth. The noise r\ is uncorre- 
cted Gaussian white noise in both space and time with 
zero average i.e., (r)(x,t)) = and (r](x,t)r](x' ,t')) = 
2DS d (x - x')S(t - t'). The KPZ equation is invariant 
under translations along both growth direction and per- 
pendicular to it, as well as time translation and rotation 
[131 ]. Rescaling the variables, h = h\j2D jv, and t = i/v, 
changes Eq. (J) to, dh(x,t)/di = V 2 h + yfe \ Vh | 2 
+fj{x,t), where e = \ 2 D/2v 3 and (fj(x,t)fj(x' ,t')) = 
S d (x — x')S(t — P). In the following, we work with the 
single parameter e and drop all the tildes for simplicity. 

We have studied the rescaled KPZ equation on a square 
lattice with periodic boundary conditions. The numeri- 
cal integration was done using the Runge-Kutta-Fehlberg 
scheme of orders 0(4) and 0(5) [l4[. This scheme con- 
trols automatically the integration time step 5t, such 
that the resulting height error 5h (which is estimated 
by comparing the results obtained from the 0(4) and 
0(5) integrations) can be ignored at each time step. We 
took the error to be less than 0.1, and we checked that 
smaller values of Sh do not improve the precision of the 
computed quantities. The noise r\ was generated by the 
Box-Muller method. To avoid the instabilities that may 
appear during the growth, we used the algorithm intro- 
duced in [IH . where the term (1 — c~ 1 e~ c f), is used in- 
stead of the nonlinear term in Eq. ((T|), i.e., / =| V/i | 2 . 
Since / <C u>|(oo), where u>z,(t) is the interface width of 
the system with size L at time t, by keeping c <C w£ 2 (oo), 
one can control the possible numerical divergencies that 
may appear during the integration. Clearly for very small 
c this term converges to /. 

We have checked that the growth exponents for (1 + 1)- 
dimension are obtained correctly (both roughness and 
growth exponents a = 1/2 and (3 — 1/3 respectively), 

and the (2+l)-dimensional results are given in Fig^ 1, 

which are in good agreement with previous studies [161 ] . 

We now consider the saturated 2D-KPZ surface and set 
its mean height to be zero, and attribute the same sign 
to the points that have positive or negative heights. The 
same-sign regions (clusters) and their boundaries (loops) 
were identified by the Hoshen-Kopclman algorithm (Fig. 
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FIG. 1: (color online). Main frame: Interface width wi,(t) vs 
time t of the KPZ equation in (2+l)-dimensions and for e = 
10, for different square lattice size L. The slope of the straight 
line yields the growth exponent j3 = 0.23 ± 0.01. Upper-left 
inset: Saturation width u)l(oo) for systems of different size L. 
The slope of the solid-line fit yields the roughness exponent 
a — 0.37 ± 0.01. Lower-right inset: Rescaled wl vs rescaled 
t. 




FIG. 2: (color online). The clusters with positive heights 
are shown for the 2D-KPZ interface with different colors. 
Negative-height regions are colored with black. 



2). 

To investigate the scaling behavior of such loop ensem- 
bles in the 2D-KPZ interface, a set of scaling exponents 
associated with cluster and loop statistics were computed 
and checked, and are shown in Fig. [3] [H GJ, [ll • The 
estimated exponents are in good agreement with the cor- 
responding analytical results for SAWs in the scaling 
limit. The fractal dimension Df of contour lines ob- 
tained from the scaling relation between their length I 



and radius of gyration R, i.e., I ~ R D s ; is given by 
Df = 1.33 ± 0.01 which is in agreement with the one 
obtained by the box-counting method for the largest con- 
tour lines, Df = 1.33 ±0.02. Comparing with the known 
fractal dimension of SLE K curves Df = 1 + k/8 for 
< k < 8 , the contour lines may have conformal in- 
variant scaling limit according to SLE s / 3 . 

The quantity that can confirm the SAW property of the 
contour lines is the restriction property. Suppose that S 
is a hull in the upper-half plane H which is bounded away 
from the origin, and 7 is a simple SLE curve in H with 
k < 4. Let \&s be a unique conformal map of M\S onto H, 
such that * s (0) = 0, * s (oo) = 00 and *s(oo) = 1. The 
restriction property states that the distribution of curves 
conditioned not to hit S is the same as the distribution 
of curves in the domain M\S. This happens only for 
k = 8/3, and it is shown that @ the probability that a 
curve does not hit the hull S is 

p[ 7 ns = 0] = |*s(0)| 5/s . (2) 

To examine this property directly for contour lines on the 
saturated 2D-KPZ surface, we proceed as follows. First, 
we identify all the cluster boundaries (contour lines): for 
each cluster an explorer walks on the zero height line as 
keeping the sites with positive height on the right. Then, 
we consider an arbitrarily placed straight line for each 
curve as a real axis and cut the portion of the curve 
above it. Using this procedure, we obtain an ensemble of 
contour lines in the half-plane which start at the origin 
and end on the real To obtain curves whose 

size is of order one, we rescale them by a factor of N", 
where v = 1/Df. Second, we consider the hull S as 
a slit placed at various distances £ from the origin and 
various heights h, for which the map *f?s is defined by, 
^s(z) = £ + y/( z — £) 2 + After mapping the curves 
by ip(z) = x 00 z/(x 00 — z) [23|, we have checked Eq. ([2]) 
for the contour lines of the 2D-KPZ surface. As shown in 
Fig. [31 the result is consistent with Eq. ((2]) and implies 
the connection between the contour lines and both the 
SAWs and SLE s/3 . 

Since the restriction property only holds for SLE 8 / 3 
curves, we test the probability that an SLE curve passes 
to the left of a given point z = pe l9 , where 9 is the angle 
between the point and the origin, and p is the distance 
from the origin inside the upper half-plane. Given scale 
invariance, this probability is independent of p, and the 
theory of SLE predicts [2(| that 

p ^ =\ + v^m 2Fi G- ^ ; - cot2{6) ) com 

(3) 

Here, 2-F1 is the hypergeometric function. The computed 
P' K {0) for the contour lines of the 2D-KPZ surface is also 
consistent with the analytical form with k = 8/3± 1/10 

(Fig. HD. 

These results strongly suggest that the iso-height lines 
might be, in the scaling limit, conformally invariant, giv- 
ing rise to the SLE curves with k = 8/3. To examine this 
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FIG. 4: (color online). Main frame: The probability that 
a contour line of 2D-KPZ surface in the upper half-plane 
does not hit the slits of height 0.05 < h < 0.5 placed at £ 



from origin on the real axis vs 1^^(0)1 = \£,\(£, 2 + h 2 



-l/2_ 

The solid line is the corresponding analytical prediction for 
SLEg/ 3 curves. The error bars are almost the same size as 
the symbols in the scaling region. Inset: The probability that 
such a contour line passes to the left of a point z — pe l6 , for 
p = 0.05, 0.1, 0.15, 0.2, 0.25 and 0.3. The solid line shows the 
prediction of SLE for n = 8/3. 
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FIG. 3: (color online). Cluster and loop statistics for the iso- 
height lines of 2D-KPZ surface and the EW surface are shown 
by the lower (different colors) and upper (green) graphs, re- 
spectively. The results for the EW have been shifted by a 
constant 1 in order to distinguish them from the 2D-KPZ re- 
sults, (a) The average area M vs the radius of gyration R. 

(b) The length of a loop I vs the radius of gyration R. (c) 
Number of clusters of area between M and AM. (d) Number 
of loops of length between I and XL (e) Number of loops of 
radius of gyration between R and XR. (f) Number of loops 
of area between A and XA. (g) The average area of loops A 
vs the length I. In all figures A ~ 1.05. Solid lines show the 
corresponding analytical results for the SAWs (bottom) and 
0(2) model (top). The error bars are almost the same size as 
the symbols in the scaling regions, and have not been drawn. 
The change in the exponents relative to the SAWs in Figs. 

(c) , (d), (e) and (f), is due to the roughness exponent [13]. 



suggestion directly, we can extract the Loewner driving 
function £ of the curves using the successive conformal 
maps. We use the algorithm introduced by Bernard et al. 
based on the approximation that driving function is a 
piecewise constant function. Each curve is parameterized 
by a dimcnsionlcss parameter t, to be distinguished from 
time in {!]) . The procedure is based on applying the map 



G t ,c = x O0 {'qx 00 {x 00 -z) + [x i 00 {z-Tii) 2 +At( : x 00 -z) 2 {x 0C - 

ri?] 1/2 }/{xUx°o -*) + ~ V) 2 + ±t( Xoo - z) 2 { Xao - 

T?) 2 ] 1 / 2 } on all the points z of the curve approximated by 
a sequence of {zq — 0, z±, ■ ■ ■, Zjv = 2:00} in the complex 
plane, where r\ = </3 -1 (C) and again ip(z) — XooZ / (xoo—z). 
At each step, by using the parameters 770 = ^^(Co) = 
\Rez\Xoo — (Rezi) 2 — (Imzi) 2 ]/(x oc — Rez\) and t\ = 
(Imzi^x^/iiKRezi - Xoo) 2 + (Imzi) 2 } 2 }, one point of 
the curve zq is swallowed and the resulting curve is rear- 
ranged by one element shorter. This operation yields a 
set containing numbers of Cfc(^fc) f° r each curve. The 
next step is analyzing the ensemble of the driving func- 
tions C(t) which can indicate, within the statistical errors, 
whether the curves are SLE or not. As shown in Fig. 
the statistics of the ensemble of converges to a Gaus- 
sian process with variance (( 2 (t)) = nt and k = 2.6 ±0.1. 
This evidence certifies that the iso-height lines of 2D- 
KPZ interface in the saturation regime appear to be con- 
formally invariant and are described by the SLE s / 3 . The 
above results were obtained for e = 10; however, the same 
analysis for growth surfaces with other values of e (which 
were checked for e = 5 and 25) indicates no changes. 

In the case of e = 0, which corresponds to the EW 
model, comparing Fig. [6] and Fig. [2l indicates more 
"porosity" in the clusters, which is indicative of changes 
in the cluster boundaries' shape. As presented in Fig. 
[31 the cluster and loop statistics in this case are most 
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FIG. 5: (color online). Statistics of the driving function C(t). 
Main frame: the linear behaviour of (C(i) 2 ) with the slope 
K = 8/3 ± 1/10. Lower-right inset: The probability density 
of rescaled by its variance k at times t — 0.01, 0.015, 0.02 
which is Gaussian. Upper-left inset: The correlation function 
of the increments of the driving function S((t). Kolmogorov- 
Smirnov (K-S) goodness of fit for normal distribution of the 
noise ((t)/VKt for k = 8/3 and < t < 0.05 yields K -S = 
0.017 



consistent with those for the 0(2) model. These lead 
to the conclusion that if one assumes that the scaling 
limit of such contour lines exists, it should belong to the 
SLEi curves. We also checked this directly as above and 




FIG. 6: (color online). Positive height connected domains of 
2D-KPZ interface without the non-linear term, corresponding 
to the EW model. Negative height regions are black. 



found that the driving function has Gaussian statistics 
with variance n = 3.7 ± 0.2. 

The height clusters and their boundary statistics, in 
the manner presented here, can be applied to model the 
experimental grown surfaces by using an ensemble of the 
grown samples in the saturated regime. Since this analy- 
sis is far more accurate, it may be also used to investigate 
wether a model belongs to a universality class or not. For 
example, the small difference between the cluster analy- 
sis of the BD model @, HO, S3 and KPZ equation in two 
dimensions can be revealed. 
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work. We also thank S. Moghimi-Araghi, M. A. Ra- 
jabpour and M. Sahimi for useful discussions. 
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